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I.  OBJECTIVE 

The  work  described  in  this  report  is  aimed  at  demonstrating  the 
feasibility  of  the  Nuclear  Thermionic  Radiator  Space  Power  System.  It 
is  a  continuation  of  work  initiated  on  Contract  AF  33(616) -8 119 .  The 
specific  objectives  of  the  current  program  are: 

1)  To  demonstrate  operation  of  at  least  three  thermionic 
converters  in  series  array,  using  heat  from  a  liquid  metal 
loop  to  heat  the  cathodes 

2)  To  demonstrate  1000  hours  of  thermionic  operation  using 
liquid  metal  heating 


o 


l 


AN- 97 3 


II.  INTRODUCTION 

The  thermionic  radiator  system  will  utilize  a  nuclear  heat  source  to 
produce  electrical  power  in  space.  A  combined  thermionic  converter-radiator 
unit  is  coupled  with  a  high  temperature,  liquid  metal  cooled,  fast  reactor 
heat  source;  coolant  is  circulated  directly  through  the  radiator  and  back 
to  the  reactor  in  a  single  loop.  Cylindrical  thermionic  converters  filled 
with  cesium  vapor  encircle  the  tube  walls  of  the  radiator.  Emitters  (on  the 
radiator  tubes)  are  heated  by  conduction  through  the  Lube  wall.  Collect¬ 
ors  concentric  to  the  cathode  dissipate  waste  heat  directly  to  space 
through  attached  fins. 

One  aspect  of  the  feasibility  of  the  thermionic  radiator  system  is  to 
be  demonstrated  by  operating  three  generators  with  heat  from  a  single  liquid 
metal  loop.  The  generators,  to  be  electrically  connected  in  series,  are  of 
a  type  being  developed  by  RCA  under  Contract  AF  33(657) -8005 .  The  con¬ 
verters  are  cylindrical  and  will  be  mounted  on  a  columbium  alloy  tube,  using 
an  electrically  insulating  thermal  bond.  The  liquid  metal  loop  on  which 
these  converters  will  be  mounted  will  be  fabricated  from  the  same  columbium 
alloy.  Lithium  will  be  circulated  in  the  loop  by  natural  convection.  The 
entire  system  will  be  placed  in  a  vacuum  chamber  to  protect  the  refractory 
metals  from  oxidation. 

In  addition  to  this  experimental  work,  a  thermal  transport  analysis 
is  being  conducted  to  aid  in  selection  of  loop  test  parameters  and  to  allow 
overall  system  characteristics  to  be  investigated.  The  current  status  of 
the  work  is  as  follows: 
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1)  The  thermal  transport  analysis,  except  for  details  of  radiator 
configuration,  has  been  completed. 

2)  The  liquid  metal  loop  is  partly  constructed.  Delivery  of  the 
thermionic  generators  has  been  delayed  because  of  fabrication 
difficulties  associated  with  the  electrically  insulating  thermal 
bond  between  the  generator  and  the  liquid  metal  tubing. 
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III.  TECHNICAL  PROGRESS 

Technical  work  on  this  program  has  been  divided  into  the  following 
three  tasks: 

Task  2  -  Liquid  Metal  Tubing  -  Converter  Bond 

Task  3  -  Thermal  Transport  Analysis 

Task  4  -  Series  Unit  Loop  Test 

Task  2  was  completed  prior  to  the  period  covered  in  this  report.  Technical 
progress  on  Tasks  3  and  4  is  discussed  below. 

A.  TASK  3  -  THERMAL  TRANSPORT  ANALYSIS 

During  this  quarter,  attention  has  been  directed  to  degradation 
of  converter  performance  resulting  from  series  operation,  the  effect  of 
emitter  temperature  variation  on  converter  performance,  and  calculation  of 
micrometeoroid  armor  requirements.  These  separate  efforts  are  described 
below . 

1 .  Effect  of  Emitter  Temperature  on  Converter  Performance 

The  effect  on  power  density  of  series-connecting  converters 
operating  at  different  emitter  temperatures  for  the  range  1300  to  1175°C  was 
found  to  result  from  connecting  converters  electrically  in  series  over  this 
temperature  region.  Non-optimum  performance  can  result  from  connecting 
converters  electrically  in  parallel  when  the  emitter  temperatures  differ, 
unless  output  voltages  are  carefully  matched.  The  effect  on  power  density 
of  a  linear  emitter  temperature  variation  in  an  individual  converter  was 
calculated  for  a  temperature  drop  of  up  to  40°C. 
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Power  and  efficiency  can  be  lost  when  thermionic  converters 
are  series-connected,  if  units  are  not  perfectly  matched  with  respect  to  the 
current  produced  at  optimum  operating  conditions.  Mismatching  can  occur 
either  from  variation  in  fabrication  and  processing  of  the  particular  con¬ 
verters,  or  from  variation  in  operating  temperature;  only  this  second 
consideration  is  discussed  here.  Since,  in  the  thermionic  radiator  approach, 
heat  is  extracted  from  a  liquid  metal  stream  to  supply  converters  with 
thermal  energy,  the  cathode  temperature  will  decrease  along  the  length  of  a 
series  of  converters  as  well  as  along  each  individual  converter,  and  some 
degradation  in  output  will  result  from  this  temperature  differential  between 
series-connected  units. 

a.  Series  Connection 

To  evaluate  the  effect  of  series-connecting  thermi¬ 
onic  converters,  representative  current-voltage  performance  characteristics 
must  be  obtained.  The  shortage  of  experimental  data  at  the  low  temperature 
operating  range  of  interest  (about  1200°C)  severely  handicaps  this  analysis. 
Although  it  was  intended  to  collect  and  compare  all  published  data  for  the 
molybdenum  cathode-anode  system  to  determine  characteristic  voltage-current 
performance,  lack  of  reproducibility  obtained  by  various  investigators 
testing  converters  under  operating  conditions  stated  to  be  identical,  made 
it  impossible  to  establish  such  information;  it  was,  therefore,  necessary 
to  use  data  from  a  single  source.  Calculations  are  based  on  the  findings 
of  Baum  and  Jensen  (General  Electric  Company)  shown  in  Figure  1.  The 
results  of  extrapolating  this  information  to  lower  temperatures  over  a 
limited  I-V  region  are  given  in  Figure  2.  Low  voltage  arc  mode  operating 
characteristics  for  a  range  in  emitter  temperature  of  1175  to  1300°C  are 
summarized  in  Figure  3. 

The  series-connection  of  converters  demands  that  the 

current  be  the  same  in  each  device.  As  evident  from  Figure  3,  this  con- 
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dition  can  be  satisfied  (about  7  amps/cm  )  throughout  the  specified  temper¬ 
ature  range  without  deviating  significantly  from  the  maximum  power  condition 
for  each  operating  temperature;  fortuitously,  peak  power  is  obtained  at 
almost  the  same  current  density.  At  high  cesium  pressure,  variations  in 
converter  performance  attributable  to  temperature  differences  are  minimized 
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by  a  self-compensating  effect.  A  temperature  decrease  lowers  the  emitter 
work  function  due  to  increased  cesium  coverage;  this  tends  to  increase 
electron  emission  and  compensate  for  the  loss  caused  by  lower  emitter 
temperature.  This  compensation  is  very  nearly  complete  for  the  case  con¬ 
sidered.  As  a  result  of  the  lowered  emitter  work  function,  less  output 
voltage  is  obtained,  and  the  net  effect  of  a  temperature  decrease  is  less 
power. 

For  the  operating  region  considered,  no  significant 
degradation  in  converter  performance  due  to  non-optimum  operation  is  ex¬ 
pected  as  a  result  of  series-connecting  converters  operating  at  different 
temperatures.  Near  peak  performance  can  be  obtained  from  each  converter; 
however,  the  power  output  from  each  device  will  obviously  decrease  as  the 
emitter  temperature  decreases.  The  effect  of  emitter  temperature  on  peak 
power  density  is  shown  in  Figure  4.  Because  converter  power  density  is 
very  nearly  linear  over  this  temperature  range,  the  average  power  density 
of  a  number  of  series  connected  devices  corresponds  to  that  of  the  average 
emitter  operating  temperature,  since  the  deviation  from  optimum  conditions 
at  each  converter  is  negligible. 

b.  Parallel  Connection 

In  the  preceeding  analysis,  it  was  found  that  near 
optimum  performance  can  be  obtained  from  converters  connected  in  series, 
since  peak  power  is  obtained  at  almost  the  same  converter  current  density 
over  the  operating  temperature  range  considered  (1175  to  1300°C) .  This  is 
apparent  from  the  converter  operating  characteristics  shown  in  Figure  3. 

It  is  equally  apparent  that  sizeable  divergence  from  optimum  conditions  can 
be  incurred  by  operating  converters  in  parallel  where  the  voltage  across 
each  device  is  the  same.  For  example,  if  the  voltage  corresponds  to  maximum 
power  for  an  emitter  temperature  of  1300°C  (0.29  volts) ,  a  converter  con¬ 
nected  in  parallel  operating  with  an  emitter  temperature  of  1175°C  would 
provide  negligible  power  output.  This  serious  consequence  is  emphasized 
by  considering  the  parallel  connection  of  two  individual  converters 
operating  at  significantly  different  temperatures.  Power  loss  from  con¬ 
verters  operating  in  parallel  can  be  completely  eliminated  by  series¬ 
connecting  a  sufficient  number  of  low  temperature  converters  to  match  the 
optimum  voltage  output  of  the  high  temperature  converters.  In  the  thermionic 
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radiator  system,  a  number  of  series-connected  converters  are  joined  in 
parallel  to  similar  radiator  segments,  each  operating  over  the  same  total 
temperature  difference.  Since  temperatures  are  identical,  no  power  loss 
will  result  from  parallel  connections,  if  each  converter  can  be  fabricated 
to  provide  identical  operating  characteristics.  But  if  large  nwnbers  of 
series  connected  converters  show  significant  variation,  it  may  be  necessary 
to  increase  flexibility  during  assembly  so  that  the  proper  number  of  con¬ 
verters  can  be  series-connected  in  each  radiator  segment  to  assure  matched 
total  output  voltage. 

c.  Temperature  Variation  in  Individual  Converter 

Thus  far,  the  analysis  is  based  on  an  individual 
converter  operating  at  a  uniform  emitter  temperature;  however,  since  heat 
is  continuously  removed  from  the  liquid  metal  stream,  the  emitter  temper¬ 
ature  will  vary  along  the  length  of  each  converter.  The  observed  converter 
output  should  be  equivalent  to  that  of  a  large  number  of  diodes  connected 
electrically  in  parallel,  but  operating  at  different  emitter  temperatures, 
as  illustrated  in  Figure  5. 


EMITTER 


COLLECTOR 


M<  T1 


FIGURE  5.  EQUIVALENT  CIRCUIT  FOR  CONVERTER  WITH 
EMITTER  TEMPERATURE  VARIATION 

From  the  preceding  considerations,  it  is  apparent  that,  if  the  temperature 
variation  is  large,  low  temperature  emitter  regions  of  the  converter  may 
be  very  ineffective  and  that,  consequently,  the  temperature  difference  along 
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an  individual  converter  must  be  minimized.  Since  the  total  temperature 
drop  along  a  single  converter  will  depend  on  the  overall  radiator  configu¬ 
ration,  but  in  no  case  will  exceed  40°C ,  the  power  loss  associated  with  an 
emitter  temperature  variation  can  be  estimated  as  follows.  A  typical 
motive  diagram  for  high  cesium  pressure  thermionic  converter  operation 
(Figure  6)  will  be  used  as  a  model  for  the  analysis. 


FIGURE  6.  TYPICAL  THERMIONIC  CONVERTER  MOTIVE  DIAGRAM 
(HIGH  CESIUM  PRESSURE  OPERATION) 


The  net  current  flowing  through  the  converter  can 
be  expressed  at  the  emitter  side  in  terms  of  the  potential  V„  of  the  emitter 
sheath  and  the  random  electron  current,  1^,  in  the  plasma: 


I  =  I 

s 


I 

e 


e 


-e  V^/k  T 
E  e 


(1) 


The  same  current  can  be  expressed  at  the  col  lector  side  in  terms  of  the 
potential  of  the  collector  sheath  and  the  random  electron  current  in  the 
plasma :* 


*E.  N.  Carabateas,  S.  D.  Pezaris,  and  G.  N.  Hatsopoulos,  "Intepretation  of 
Experimental  Characteristics  of  Cesium  Thermionic  Converters,"  J ■  Appl  ■ 
Phys.  ,  32  ,  352  (1961). 
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where 
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I 

e 

I 
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c 

T 


e 


I 

e 


-e  V  /k  T 
e  c  e 


net  current, 

random  electron  current, 
emitter  saturation  current , 
emitter  sheath  potential, 
collector  sheath  potential,  and 
temperature  of  electrons. 


(2) 


The  assumption  was  made  that  the  plasma  is  uniform,  which  implies  that  the 
electron  concentration,  electron  temperature  and,  therefore,  random  electron 
current  are  the  same  at  the  two  edges  of  the  plasma  next  to  the  emitter  and 
collector.  The  ion  current  was  neglected  on  both  the  emitter  and  collector 
side,  since  it  is  generally  small  compared  to  the  net  electron  current. 

From  Equation  (2)  , 


I 

e 


Ie 


V  /k  T 
c  e 


and  by  eliminating  in  Equation  (1) , 


I  =  I  -  Ie 


<VE- 


V  )/k  T 
c  e 


solving  for  I,  then, 


I  -  I  /!  +  .-'  <#E  -  V/k  Te 

s 


(3) 


Assuming  that  the  collector  and  emitter  sheath  potentials  and  electron 
temperature  are  constant  for  constant  output  voltage  over  the  small  emitter 
temperature  range  considered, 


1  +  e 


-e  (V 


V  )/k  T 
c  e 


*  e 


constant , 


the  net  electron  current  density,  J,  can  be  expressed  as: 

J  =  AT2  e‘e  (VK)/kI  , 


(*) 
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where  A  =  Richardson  constant, 

T  =  emitter  temperature,  and 

0  =  emitter  work  function. 

E 

The  observed  total  output  current  from  the  converter 
is  then  the  sum  of  the  currents  from  each  of  the  small  diodes  connected  in 
parallel  along  the  emitter  surface.  The  average  current  density,  J,  for  a 
cylindrical  emitter  can  therefore  be  expressed  as 


J(X)  dx,  (5) 


where  J(x)  =  A  T  (X)  exp 

and  0  =  0  +  K  =  constant. 

E 

The  emitter  work  function,  0  ,  is  assumed  to  be 
constant  over  the  small  temperature  range  considered,  which  is  only  strictly 
valid  for  vacuum  converters.  This  represents  a  worst  case  solution,  then, 
inasmuch  as  the  emitter  work  function  decreases  as  the  temperature  decreases, 
and  would  tend  to  reduce  the  calculated  emission  loss.  For  a  linear  emitter 
temperature  distribution, 

T (X)  =  Tq  (1  -  OCX)  , 

the  integral  equation  can  be  expressed  as: 

at2  rro 
3  =— 

«  L  J  0 
1 


e  0 
k  T  (X) 
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where 


if  -  1/1  -  ax, 

K  •  VTL  • 


Tq  =  maximum  emitter  temperature, 

T,  =  minimum  emitter  temperature,  and 
L 


L  =  converter  length. 
The  solution  to  Equation  (6)  is: 


A  T 


J  = 


<To  '  V 
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S’ 

e  <6\ 

tl1 

3 

e  0  1 

E4 

k  T 

T 

E4 

k  T 

0/ 

oJ 

ll 

(7) 


r<*>- t 

3  G 

where  E^(X)  =  X  — dt 


Equation  (7)  was  solved  for  T  =  1200°C  and  0  <  T  -  Tt|  <  40°C ,  assuming 

0  \  O 


J  (1200°C)  =  7.0  amps/cm^; 


7  2 

therefore,  0  +  K  =  0  =  2.21  for  A  =  120  amp/cm  °K  . 

£ 

These  results  are  illustrated  in  Figure  7.  Since  the  output  voltage  is 
constant  across  the  converter,  Figure  7  also  depicts  the  effect  of  temper¬ 
ature  drop  on  converter  power  density. 

2 .  Micrometeoroid  Shielding  Calculations 

Problems  associated  with  the  micrometeoroid  shielding  of 
the  thermionic  radiator  system  were  briefly  reviewed,  including  the  latest 
estimates  of  the  micrometeoroid  flux  distribution  and  comparisons  of  various 
penetration  models.  Equations  describing  the  required  armor  thicknesses 
were  formulated.  The  review  and  derivations  appear  as  Appendix  1.  The 
results  of  this  analysis  can  be  expressed  in  terms  of  two  equations  for 
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To  =  1 200°C 
J  (1200)  =  7  amp/cm^ 


To  -  Tl  (°C) 


FIGURE  7.  EFFECT  OF  LINEAR  EMITTER  TEMPERATURE  DROP 
ON  CONVERTER  CURRENT  DENSITY 
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the  equivalent  armor  thickness.  The  protection  afforded  the  liquid  metal 
tubing  by  converter  materials  and  armor  is  given  by: 


t 

equiv 


/_  \ 0. 334  /  \ 

E 

tube 

I  ^tube  \ 

E 

fi 

armor 

1  ^armor  / 

0.166 


Ctube 


where 


E  \0'334 

/o 

conv  1 

conv 

E 

fi 

armor/ 

armor 

t 

conv 


+  t 


external  armor 


E  =  Young's  Modulus  of  material  indicated, 
fi  =  density  of  material  indicated,  and 
t  =  thickness  of  material  indicated. 


The  armor  thickness  required  for  a  given  radiator  area,  assuming  a  specified 
mission  time  and  the  probability  that  no  punctures  will  occur,  is  given  by: 


0.249  0.249 


5.5 

r 

1.  ,  '.Ofl  te<luiv| 

„  0.334  „  0.166 

E  /O 

L  a  a  J 

-  In  P(o) 

N  [Ai  +  2^L^48| 

where 


Ea  =  Youngs  Modulus  of  armor, 

=  density  of  armor, 

T  =  mission  time, 

P(o)  =  probability  of  no  critical  damage, 

A^  =  inside  surface  area  of  radiator  tube 
L  =  length  of  radiator  tube,  and 
N  =  Number  of  radiator  tubes. 


Using  these  relations,  armor  calculations  were  made  for  the  following 
parameters . 
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Armor  Parameters 
Material 
Density 

Young's  Modulus 

Tube  Parameters 
Material 
Density 

Young's  Modulus 
Wall  Thickness 

Converter  Parameters 
Material 
Density 

Young's  Modulus 

Wall  Thickness  (cathode  +  anode) 

Penetration  Parameters 
Mission  Time 

Probability  of  Non-puncture 


Beryllium 
1-85  gm/cc 
37  x  106  psi 


Cb-lZr 
8 .56  gm/cc 
14.4  x  10  psi 
0.030  in.  =  0.076  cm 


Molybdenum 
10.2  gm/cc 
42  x  10^  psi 
0.040  in.  =  0.102  cm 

104  hrs  =  417  days 
957. 


Substitution  of  these  parameters  into  Equations  (1)  and  (2)  yields  the 
following  results: 


t 

equiv 


+ 


14.4! 

10.334 

[8.5 

37  J 

1 

1.8 

N° 

.334 

f  10.2 

l37) 

[l.85 

0.166 


0.076  + 


0.166 


0.102  +  t 


=  0.0715  +  0.141  +  t  =  0.213  +  t 
a  a 


(3) 


equiv 


5.5 


(37xl06)0,334  ( 1 .85) °* 166 


417 


-|0.249 


L-  In  0.95J 


N  (A.  +  0.206  L  t  ) 
l  eq7 


0.249 


t 

eq 


0.01485 


9.4 


1  0.249 


N(A.  +  0.206  L 


teq^ 
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=  0.1395 


N  (Af  +  0.206  L 


0.249 


Combining  Equations  (3)  and  (4)  and  simplifying  gives: 


(4) 


(NA.)  +  0.206  (NL.)  S  =  2730  6  4-02  ,  (5) 

<3  *  0  213 

where  t  (in)  =  - - .  If  tg  is  negative,  then  no  armor  is  required 

beyond  that  provided  by  the  tube  and  converter.  The  solution  for 
Equation  (5)  for  a  selected  range  of  variables  is  shown  in  Figure  8. 

Table  1  provides  an  example  of  the  application  of  this 
curbe  for  the  following  conditions: 

Power  level  =  300  kw(e)  net  output 

=  345  kw(e)  output  at  radiator 

Emitter  diameter  =0.8  inch 

Tube  length  required  for  seals  =  5%  emitter  length 


TABLE  1 

MICROMETEOROID  SHIELDING  REQUIREMENTS 


2 

Emmitter  power  density  -  w/cm 

1.0 

0.5 

2.0 

2 

Total  emitter  area  -  ft 

372 

248 

186 

2 

Total  wetted  area  -  ft 

342 

228 

171 

Total  tube  length  -  ft 

1870 

1250 

933 

Shielding  required  -  inch 

0.19 

0.16 

0.14 
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B.  TASK  4  -  SERIES  UNIT  LOOP  TEST 

The  natural  convection  loop  described  in  the  previous  quarterly 
report  has  been  partly  fabricated.  The  general  loop  configuration  is  shown 
in  Figure  9,  and  the  construction  status  of  major  loop  components  is 
described  below: 

1.  Main  Heater 

All  molybdenum  and  columbium  heater  components  have  been 
fabricated.  The  main  molybdenum  bar  was  machined  from  sintered  molybdenum; 
the  studs  and  nuts  were  machined  from  arc-cast  molybdenum.  The  weight  of 
the  molybdenum  bar  will  be  transmitted  to  the  support  stand  by  a  spring- 
loaded  support  rod,  and  spring  tension  will  be  adjusted  so  that  the  columbium 
loop  will  be  unloaded  at  operating  temperature. 

2 .  Heat  Sink 

Both  heat  sink  units  have  been  completely  fabricated. 

All  braze  joints,  except  those  between  the  copper  caps  and  the  stainless 
steel  tubes,  were  made  with  vacuum  brazing  techniques.  Braze  alloys  and 
temperatures  are  summarized  in  Table  2.  The  stainless  steel  tubes  will  be 
brazed  to  the  copper  heads  by  means  of  a  tungsten/ inert  gas  brazing 
technique  applied  in  an  inert  atmosphere  chamber.  Successful  test  brazes 
have  been  made  on  a  heat  sink  mockup,  using  a  silver-copper-lithium  alloy. 


TABLE  2 

HEAT  SINK  BRAZE  SUMMARY 


Joining  Materials 

Columbium-Molybdenum 

Molybdenum-Copper 

Copper-Copper 


Braze  Alloy 

Gold-Palladium 
Gold-Nickel 
Si lver-Copper 


Braze  Temperature 

2200°F 

1750°F 

1430°F 
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3.  Support  Stand 

The  stainless  steel  support  stand  on  which  the  loop  will 
be  mounted  has  been  fabricated.  The  brackets  which  will  support  the  loop 
at  various  points  have  also  been  constructed  and  include: 

1)  A  rigid  support  for  the  upper  heat  sink 

2)  A  vertical  support  for  the  lower  heat  sink 

3)  The  spring-loaded  support  for  the  molybdenum 

heater  bar 

4)  The  thermionic  assembly  support 

This  last  support  is  attached  to  the  loop  at  the  top  of  the  thermionic 
assembly.  The  support  beam  is  counter- balanced  to  compensate  for  the 
weight  of  the  three  thermionic  generators,  thus  eliminating  the  bending 
load  on  the  section  of  the  loop  between  the  thermionic  generators  and  the 
heat  sink. 


4.  Thermionic  Generators 

The  thermionic  generators  for  this  loop  are  being 
fabricated  by  RCA  under  subcontract  to  AGN ,  and  are  identical  to  those 
developed  by  RCA  for  ASD  under  Contract  AF  33(657) -8005.  The  general 
specifications  of  the  converter  are: 

1)  The  converter  emitter  will  be  molybdenum. 

2)  The  center  mounting  tubing  will  be  columbium- 17, 
zirconium. 

3)  The  entire  converter  assembly  will  be  capable  of 
withstanding  a  temperature  of  500°C  without  deleterious 
effect . 


4)  The  converter  will  employ  a  concentric  cylindrical 
configuration,  as  follows: 


Emitter  0D 

Emitter  active  length 

Emitter  area 

Converter  length 
(excluding  tubing) 


0.85  to  1.050  inches 

2.0  to  2.5  inches 
2 

40  cm 

4.25  inches  (maximum) 


I 

i  i 
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Figure  10  is  a  drawing  of  the  complete  converter  assembly. 

Four  converters  under  construction  at  RCA  were  scheduled 
for  delivery  to  AGN  by  11  April  1963.  Just  prior  to  shipment,  however,  a 
vacuum  leak-check  revealed  small  leaks  in  the  inner  columbium  tube. 
Although  the  cracks  in  no  way  affected  converter  performance,  they  would 
permit  liquid  metal  to  escape,  thus  rendering  them  useless  for  their 
intended  application.  Fabrication  was  suspended  and  the  design  is  being 
re-examined . 
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IV.  FUTURE  PLANS 

Delivery  of  the  modified  thermionic  converters  is  not  expected  until 
31  August  1963.  Loop  fabrication  and  construction  of  auxiliary  components 
will  proceed  at  a  reduced  level  of  effort  during  the  next  three  months. 

After  the  converters  have  been  delivered,  loop  fabrication  will  be  completed 
and  testing  will  begin. 
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APPENDIX  I 

MICROMETEOROID  SHIELDING  REVIEW  AND  DERIVATIONS 

I.  OBJECTIVE 

This  appendix  reviews  various  micrometeoroid  shielding  problems.  The 
most  recent  estimates  of  flux  distribution,  comparisons  of  penetration  models, 
and  equations  from  which  the  required  armor  thickness  can  be  derived  are 
presented.  A  generalized  set  of  curves,  giving  the  required  armor  thickness 
as  a  function  of  radiator  size  for  an  assumed  mission  time,  survival  probabi¬ 
lities,  and  materials  selection  in  a  typical  thermionic  radiator  system,  is 
included. 

II.  DISCUSSION 

The  design  of  a  space  power  system  necessitates  consideration  of  many 
environmental  effects.  One  of  the  most  disconcerting  conditions  confronting 
the  radiator  designer  is  that  posed  by  the  micrometeoroid  flux  in  space. 
Vulnerable  radiator  areas  must  be  shielded  to  reduce  the  probability  of 
damage  or  puncture  and  consequent  loss  of  fluid  inventory,  and  this  invariably 
increases  radiator  weight.  Micrometeoroid  problems  are  discussed  below  and 
design  procedures  for  attaining  the  desired  protection  are  derived,  using  the 
latest  flux  distributions  and  penetration  correlations. 

A.  ORIGIN  AND  COMPOSITION 

It  is  currently  believed  that  90%  of  the  meteoroids  are  of  cometary 
origin  and  that  the  remaining  10%  are  fragments  of  astroids.  The  density  of 
the  cometary  meteoroids  has  been  estimated  to  vary  from  0.05  gm/cc  for  very 
porous  particles  to  3.5  gm/cc.  The  astroidal  meteoroid  particle  densities 
are  estimated  to  vary  between  3.5  gm/cc  and  8  gm/cc.  The  most  recent  reliable 
estimates  place  the  average  meteoroid  density  at  0.6  gm/cc. 
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B.  METEOROID  FLUX  DENSITY 

Various  estimates  have  been  made  for  the  meteoroid  flux  in  the 
vicinity  of  the  earth,  the  most  recent  resulting  from  the  Harvard. 
Trailblazer  experiments  performed  in  1962.  This  information  is  plotted  in 
Figure  A-l,  together  with  comparable  data  from  other  observers.  The  data 
may  be  expressed  in  the  form  of  the  following  equation; 


where,  according  to  the  data, 

oC  =  1.01  x  io‘10, 

0  =  1.34,  and 

m  =  meteoroid  mass  in  grams. 

Observations  indicate  that  the  meteoroid  flux  is  anistropic,  that  is,  a 
function  of  direction  in  space.  It  can  be  seen  from  the  radar  observations 
given  ir.  Figure  A-2  that  the  apparent  flux  is  greatest  on  the  side  of  the 
earth  facing  the  earth's  motion,  and  that  most  of  the  particles  are  traveling 
in  the  same  direction  as  the  earth.  Figure  A-3  shows  that  most  of  the 
meteoric  particles  have  trajectories  at  small  angles  to  the  plane  of  the 
ecliptic.  Based  on  overall  anistropricity,  the  radiator  vulnerability  can 
be  reduced  by  a  controlled  orientation. 

The  velocity  limits  of  particles  observed  from  the  earth  lie 
between  11  and  73  km/sec.  The  lower  limit  is  the  minimum  theoretical 
velocity  that  a  particle  would  have  if  it  fell  from  infinity  to  the  earth 
under  the  action  of  the  earth's  gravitational  force.  The  upper  limit  is 
based  on  the  assumption  that  the  earth  runs  into  particles  which  move  in 
retrograde  orbit  around  the  sun.  Figure  A-4  gives  the  velocity  distribution 
of  shower  and  sporadic  meteors.  The  average  meteoric  velocity  is  estimated 
to  be  approximately  30  km/ sec. 
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FIGURE  A-l .  METEOROID  MASS  -  FREQUENCY 
DISTRIBUTION  (FROM  ARS  2543-62) 


ANTIHELION 


FIGURE  A-2.  POLAR  DIAGRAM  IN  PLANE  OF  EARTH'S  ORBIT  SHOWING  APPARENT  NUMBER 
OF  METEOR  RADIANTS  PER  UNIT  ANGLE  PER  SECOND  (FROM  NASA  TN-D-149 


ECLIPTIC  LATITUDE  -  degrees 


FIGURE  A-3 


.  DISTRIBUTION  OF  METEOR  FLUX  ABOUT  THE  PLANE 
OF  THE  ECLIPTIC  (FROM  NASA  TN-D-1493) 
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FIGURE  A-4.  METEORIC  VELOCITY  DISTRIBUTION  (FROM  NASA  TN-D-1493) 
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C.  HYPERVELOCITY  IMPACT 


The  hypervelocity  impact  process  is  complex  and  not  yet  well 
understood,  expecially  in  the  velocity  ranges  of  interest  to  the  radiator 
designer.  To  date,  several  different  impact  regions  have  been  established 
(Figure  A-5).  In  the  undeformed  projectile  region,  the  penetration  depth 
varies  approximately  as  the  4/3  power  of  projectile  velocity,  while  in  the 
transitional  region,  penetration  does  not  seem  to  vary  strongly  with  velocity. 
In  the  fluid  flow  region,  the  target  material  behaves  as  a  fluid.  The  per¬ 
tinent  or  meteoroid  impact  region  is  uncertain  because  no  experimental  data 
are  available.  On  the  basis  of  data  accrued  with  currently  attainable 
velocities,  many  investigators  have  obtained  results  for  penetration  of 
infinite  targets  which  can  be  expressed  by: 


=  Y  d 

I* 

0 

V 

1*1 

The  penetration  theories  derived  theoretically  or  from  limted  experimental 
data  are  compared  in  Figure  A-6.  Most  investigators  found  the  penetration 
to  vary  as  the  2/3  power  of  velocity,  although  Bjork  obtained  an  exponent 
of  1/3.  The  exponent  9  is  1/3  for  a  momentum-dependent  cratering  process 
and  2/3  for  a  kinetic  energy-dependent  process.  The  penetration  calculations 
use  1/2  for  0. 


D.  PROTECTION  METHODS 

Micrometeoroid  structural  damage  can  be  generally  categorized 
as  follows,  and  each  category  requires  a  different  degree  of  protection. 

1)  Perforation  -  results  in  loss  of  working  fluid,  cesium,  etc. 

2)  Intrusion  -  energy  of  incident  particle  transmitted  to  fluid, 
possibly  causing  damage  by  shock  wave  propogation  at  local  or  distant 
points . 

3)  Surface  Damage  -  reduces  structural  integrity,  emissive  coating 
effectiveness . 

4)  Functional  Damage  -  introduces  deformations  which  could  interfere 
with  the  components  functions. 
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FIGURE  A-5.  PENETRATION  VERSUS  DENSITY  RATIO  IMPACT 
MACHINE  NUMBER  (FROM  ARS  2543-62) 
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FIGURE  A -6.  COMPARISON  AMONG  PENETRATION 
THEORIES  (FROM  ARS  2543-62) 
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5)  Spalling  -  partial  blockage  of  flow  passages  or  short-circuiting 
of  converters  by  ejection  of  material  on  back  side  of  target. 

Several  protective  methods  against  the  meteoroid  hazard  are 
theoretically  available  to  the  radiator  designer;  direct  armoring  (massive 
or  laminated),  bumper  shields,  redundant  segmenting,  self -sealing  systems, 
and  changes  of  geometry  and  operating  temperature. 

1 .  Direction  Armoring 

In  this  approach,  a  solid  mass  of  material  physically  stops 
an  impacting  particle  before  it  can  cause  serious  damage.  This  is  the  approach 
recommended  for  thermionic  radiators  because  1)  preliminary  design  data  and 
and  direction  calcula iional  methods  exist;  and  2)  circumferential  temperature 
variations  can  be  reduced  in  planar  finned  tubes.  At  the  heat  fluxes  required 
in  thermionic  converters,  large  temperature  variations  at  the  periphery  of 
the  anode  may  present  operational  problems  which  can  be  mitigated  by  increasing 
the  conduction  path  <-hrough  the  armor. 

Laminated  armor  does  not  apply  to  thermionic  radiators 
because  of  the  increased  resistance  to  heat  transfer  imposed  by  the  laminations 
laminated  armor  usually  consists  of  a  thermal  insulator  between  solid  sheets. 

2 .  Bu'iter  Shields 

In  the  bumper  concept,  a  relatively  thin  shield  of  material 
separated  from  the  surface  to  be  protected  disintegrates  an  impacting  particle 
and  diffuses  its  energy  over  a  greater  area.  The  total  thickness  of  bumper 
and  primary  armor  required  to  stop  a  given  particle  has  been  found  to  be 
approximately  one-haV:  the  thickness  for  regular  armor,  although  considerable 
variation  has  been  observed  in  different  experiments.  Note  that,  since  the 
bumper  imposes  an  additional  resistance  to  radiation,  the  effectiveness  of  this 
approach  will  depend  upon  an  evaluation  of  the  bumper’s  heat  transfer  properties 
as  well  as  its  protective  characteristics.  Bumper  thickness  is  a  function 
of  target  and  projectile  material  properties  and  projectile  size  and  velocity 
as  well  as  separation  between  bumper  and  armor.  Although  lack  of  data  excludes 
this  method  from  preliminary  consideration,  it  should  be  reassessed  when  more 
information  becomes  available.  In  thermionic  converters,  a  cesium  gap  exists 
between  the  cathode  and  anode.  At  the  required  gap  and  anode  thickness, 


A-10 


AN-973 


experimental  evidence  indicates  that  the  bumper  effect  is  negligible;  hence, 
the  unit  may  be  treated  as  solid  armor. 

3.  Redundant  Segmenting 

This  design  consists  of  a  number  of  panels  in  parallel,  a 
specified  number  of  which  will  provide  the  required  protection.  Damaged 
panels  are  isolated  from  the  system.  If  the  probabilities  of  survival  are  high, 
and  if  a  reliable  leak-detection  and  shutoff  system  can  be  developed,  redundant 
segmenting  would  be  a  weight  saving  solution.  Further  research  is  necessary 
to  determine  whether  or  not  segmenting  will  increase  system  reliability. 

4.  Self-Sealing  Systems 

Self-sealing  demands  a  minimum  amount  of  armor,  since 
punctures  are  automatically  repaired.  However,  a  satisfactory  self-sealing 
technique  remains  to  be  developed. 

5.  Changes  in  Geometry  and  Operating  Temperatures 

Increased  operating  temperatures  might  reduce  the  vulnerable 
area,  depending  on  the  influence  of  temperature  on  efficiency.  The  desirability 
of  temperature  changes  can  only  be  determined  from  extensive  system  calculations. 
Optimization  of  individual  tube  and  overall  radiator  configuration  in  terms  of 
system  weight  would  indicate  the  armor  required  to  minimize  total  weight. 

E.  METEOROID  CALCULATIONS 

1.  Derivation  of  Armor  Equations 

The  cumulate  meteoroid  flux  can  be  expressed  as  a  function 
of  particle  mass  by  the  equation, 


F 


> 


(1) 


If  N^  is  defined  as  the  average  number  of  impacts  on  vulnerable  area.  A,  in 
time,  T",  by  a  critical  mass,  mc>  then  from  Equation  (1) 


N  =  A^F  «  A  ra  P  , 
>  >  c 


(2) 


Assuming  that  the  flux  in  space  is  isotropic  (actually  it  is  anistropic; 

to  a  randomly  tumbling  vehicle,  however,  the  intergrated  flux  appears  isotopic). 
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the  flux  emitted  per  unit  area  per  unit  solid  angle  obeys  Lambert's  cosine 
relation  and 


•l 


F  cos  X  dw 

u 


tfF, 


u 


(3) 


where  Fw 
impacts , 
by 


Is  the  Isotropic  flux  In  space.  The  differential  number  of  critical 
dN^,  occuring  per  unit  solid  angle  per  unit  projected  surface  is  given 


dN^  =  F^  Ta  cos  N  d 


rr 


'Y A  cos  A  d 


(4) 


A  relationship  between  the  critical  mass,  mc>  and  the  penetration  depth  is 
needed.  On  the  basis  of  data  obtained  in  currently  achievable  velocity  ranges, 
the  depth  of  penetration  can  probably  be  expressed  by 


P 


oO 


Y  d 


(5) 


For  a  spherical  projectile,  the  diameter  can  be  written  as 


d 


[+r  -i/3 


(6) 


The  required  armor  thickness  in  terms  of  a  thin  plate  and  spalling  factor 
"a"  and  depth  of  penetration  is 


t  -  a  • 


(7) 


The  critical  damage  appears  to  depend  on  the  impact  angle;  thus,  the  velocity 
governing  the  damage  is 

v^  =  v  (cos  A  )n  .  (8) 

Combining  Equations  (5),  (6),  (7)  and  (8),  and  solving  for  the  mass  of  the 
penetrating  particle, 
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(9) 


(10) 


Substituting  Equation  (10)  into  (4)  gives: 


dN 


> 


oi 

IT 


T A  (cos 


^)3n©/S>+  1 


(11) 


By  converting  Equation  (11)  to  spherical  coordinates  and  integrating  over 
half  space,  the  final  expression  for  the  critical  number  of  impacts  can  be 
derived: 


Because  of  the  random  nature  of  the  events,  the  meteoroid  punctures  may  be 
expected  to  follow  a  Poisson  distribution,  given  by 

(N  )n  e'N> 

P(n)  - ~ 7 -  .  (13) 

11  • 

The  probability  that  no  puncture  will  occur  is  calculated  from  Equation  (13) 
by  setting  n  ■  o; 

P(o)  -  e‘N>  .  (14) 

Solving  for  N>  from  Equation  (14)  gives 

-  -lnP(o).  (15) 
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Combining  Equations  (15)  and  (12)  and  solving  for  t  gives  the  final  expression 
for  the  armor  thickness: 


\  oLT  a  ] 

\vnB  I  7 
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1/3/2 

•  (16) 


The  nomenclature  and  constants  of  Equation  (16),  based  on  the  latest  theoretical 
and  experimental  studies,  are  listed  below. 


t 

a 

if 

Pp 

PT 

0 (. 

P 

A 

P(o) 

r 

v 

c 

where 

0 

0 


required  armor  thickness  in  cm 

finite  plate  and  spalling  factor  “  1.75 

constant  in  penetration  correlation  =  2 

mean  meteoroid  particle  density  =  0.6  gm/cc 

armor  material  density  -  gm/cc 

constant  in  flux  distribution  «  1.01  x  10  ^ 

exponent  in  flux  distribution  =  1.34 

2 

external  armor  surface  -  ft 

design  probability  for  no  critical  damage  to  radiator  system 
mission  time  -  days 

average  meteoroid  velocity  =  30  km/sec 

sonic  velocity  of  the  armor  material  -  km/sec  *  2.63  x  10 ■  ®t 
Young's  modulus  at  operating  temperature  -  lb/in2.  ^ 

exponent  in  penetration  correlation  =  1/2 
exponent  in  penetration  correlation  =  2/3 
exponent  for  velocity  impact  angle  dependence  *  1 
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Substituting  the  constants  and  simplifying  Equation  (16)  yields: 


1.75  x  2 


1/3  -1/3  (0.6.\O,5lr 

> 

T 


(0-6) 


1  .01  X  io'^Ta1 


In  P(o) 


3  x  1.34 


30  |°-667 

2 

2.63  x  10'3Ve  /  P  I 

.2  x  1  x  0.667  x  1.34  x  2. 

|3  x  1.34 

4.00 

\fE 


1.01  x  10~10  T A 


=  5.5 


2  A 

L-  In  P(o) 


In  P(o) 
10.249 


0.249 


x  500 


0.667 


x  0.853 


T 1 


0.334  o  0.166  • 


,  0.334  p 
't  V 


(17) 


A  refinement  of  the  calculations  to  include  the  fact  that  a 
variable  thickness  o£  actual  armor  material  may  be  used  in  thermionic  radiators 
because  of  initial  differences  in  component  thickness  (which  are  part  of  the 
total  armor)  will  now  be  derived.  For  m  portions  having  areas,  A^,  and 
equivalent  thicknesses,  t^,  the  joint  probability  is: 


P(°)  = 


(18) 


After  taking  the  logarithm  of  Equation  (18),  the  joint  probability  becomes 

m 


•  In  P(o)  = 


i  =  1 


(19) 


Combining  Equations  (12)  and  (19)  yields  the  final  expression  for  the  variable 
armored  case; 


-  In  P(o)  =C(.T 


2 

V  A 

r3(3 1 

fvl  39® 

3n8(3  +  2 

L  Ai 

i=l 

iKi 

i 

l‘Ji 

(20) 
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2 .  Geometry  Calculations 

Assume  that  the  radiator  can  be  divided  into  a  number  of 
modules  having  the  following  configuration. 


/  .  n  1  - - - *•/ 

f* - - -  c  c  ' 


The  vulnerable  or  outside  surface  area  can  be  expressed  as; 


1T(d  +  2  t  )  1 
T  c  c 


.|»2[rr(dM+2tM)  1M 

•  »  l  lM  +  n  ‘=>  +  2fr<2  S.  S,  +  "  lo  *.>! 


^  ^total  inside  surface  +  CM  ^^^total  length  of  manifolds) 
exposed  to  fluid 


,  (21) 

where  t  is  the  sum  of  the  tube,  converter,  and  armor  thicknesses  (where 
applicable).  That  the  tube  and  converter  supply  some  protection  from  the 
meteoroids  must  be  accounted  for  in  the  armor  calculations  be  in  the  following 
manner.  In  the  derivation  of  Equation  (17),  it  was  found  that  the  armoring 
effect  of  materials  varied  as 


+  t  (2  TTl 


total  length  of  converter  carrying  tubes) 
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ta 


_ 1 _ 

„  0.334  0.166 

Et  A 


(22) 


Since  the  various  component  materials  in  series  resist  penetration,  the 
effective  protection  of  these  materials,  in  terms  of  equivalent  thickness 
of  the  armoring  material,  may  be  expressed  as: 


t 

equiv 


K» be  \ 

0.334 

'  Aube  \ 

E 

y  "Armor j 

\  armor 

E  1 

conv 

"Aonv  | 

E 

j^.  334 

/) 

armor 

i 

^armor  J 

ttube 


t 


converter 


+  t 


external 

armor. 


(23) 


To  determine  the  required  armor  thickness  without  resorting  to  the  variable 
armor  equation  (20),  it  must  be  assumed  that  t£  =  t^;  i.e.,  that  the  bare 
manifold  and  converter  thicknesses  are  the  same.  This  approximation  is 
probably  realistic,  since  a  relatively  thick  manifold  will  be  required  for 
structural  rigidity.  Using  this  assumption,  Equation  (21)  simplifies  to: 


A  =  N 


modules 


( 


A.  ,  , 
inside 


surface 


+  t  (2TL 


tubes  +  Manifolds 


'} 


(24) 


Equations  (17) ,  (23)  ,  and  (24)  can  now  be  combined  to  yield  an  expression 
for  the  required  armor  thickness  as  a  function  of  materials,  probability  of 
non-puncture,  mission  time,  number  of  modules  and  inside  surface  area,  and 
developed  length  of  each  module: 


-i 

- 

1 

0.249 

- 

5.5 

T 

N 

„  0.334  „  0.166 

LEa  A  J 

In  P(o) 

+  27TL 


tequiv 

30.48 


0.249 


(25) 
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